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NUMERICAL APPROXIMATION OF FRACTIONAL POWERS 
OF REGULARLY ACCRETIVE OPERATORS 

ANDREA BONITO AND JOSEPH E. PASCIAK 


Abstract. We study the numerical approximation of fractional powers of 
accretive operators in this paper. Namely, if A is the accretive operator as¬ 
sociated with a regular sesquilinear form A(-, •) defined on a Hilbert space V 
contained in L^(Q), we approximate A~^ for (3 € (0,1). The fractional powers 
are defined in terms of the so-called Balakrishnan integral formula. 

Given a finite element approximation space C V, A~^ is approximated 
by A^^TTh where A^ is the operator associated with the form A(-, ■) restricted 
to Yh and is the L^(Q)-projection onto V/^. We first provide error estimates 
for {A~^ — in Sobolev norms with index in [0,1] for appropriate /. 

These results depend on elliptic regularity properties of variational solutions 
involving the form A(-,-) and are valid for the case of less than full elliptic 
regularity. We also construct and analyze an exponentially convergent sine 
quadrature approximation to the Balakrishnan integral defining A^^tt/^/. Fi¬ 
nally, the results of numerical computations illustrating the proposed method 
are given. 


1. Introduction. 

The mathematical study of integral or nonlocal operators has received much 
attention due to their wide range of applications, see for instance [H1 Ian [a [22 

EHISElllidl]. 

Let 12 be a bounded domain in d > 1, with a Lipschitz continuous boundary 
r which is the disjoint union of an open set Tn and its complement F£> = F \ Fat. 
We dehne V to be the functions in H^(n) vanishing on F^i. We then consider a 
sesquilinear form A(-, •) defined for u,v gY given by 

d d 

Y, aij{x)ui{x)vj{x) +'^{ai,o{x)ui{x)v{x) 
i,j — l i—1 

+ ao^^{x)u{x)vi{x)) + aofiix)u{x)v{x) 

Here the subscript on u and v denotes the partial derivative with respect to Xi, 
i = and v denotes the complex conjugate of v. We further assume that 

A{-,-) is coercive and bounded (see, ([7]) and (jS]) below). Such a sesquilinear form 
is called regular m- 

There is an unbounded operator A on LF‘{Vt) with domain of definition D{A) 
associated with a regular sesquilinear form (see m and Section [2] below). The 
unbounded operator associated with such a form is called a regularly accretive 
operator [ST]. For such operators, the fractional powers are well defined, typically. 
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A{u,v) := / 
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in terms of Dunford-Taylor integrals. When 0 < /3 < 1, one can also use the 
Balakrishnan formula [3 [301 El]: 


( 2 ) 


A-P = 


sin(/ 37 r) 


^ + A) ^ d^. 


In this paper, we propose a numerical method for the approximation A~P f based 
on the finite element method with an approximation space V;i C V. It is worth 
mentioning that a similar formula can be used to represent the solution of a Cauchy 
problem. This is the focus of [lOj . 

Several techniques for approximating S~P f are available when S is the operator 
associated with a Hermitian form (or symmetric and real valued on a real-valued 
functional space). The most natural technique involves approximating S~P by 
where Sh is a discretization of Sh, e.g., via the finite element method using Yh- In 
this case, f for / S Yh can be expressed in terms of the discrete eigenvector 
expansion H5] ' 


= Y1 where / = ^ 

3 3 

Here denote the eigenpairs of Sh- An alternative approach is based on 

a representation of S~P f via a “Neumann to Dirichlet” map m- The numerical 
algorithm proposed and analyzed in min] consists of a finite element method 
in one higher dimension. It takes advantage of the rapid decay of the solution in 
the additional direction enabling truncation to a bounded domain of modest size. 
A third approach which is valid for more general A, is based on finite element 
approximation with an analysis employing the Dunford-Taylor characterization of 
A~Pf [231 (see, also, [311 SSj) and is most closely related to the approach which 
we will take in this paper. However, the analysis of |23j only provides errors in 
L^(fl), requires full elliptic regularity and fails to elucidate the relation between 
the convergence rate and the smoothness of /. For example, the result in [23] does 
not hold for problems on non-convex domains or problems with jumps in the higher 
order coefficients. 

The approach that we shall take in this paper is based on ([2]). The introduction 
of finite elements on a subspace Yh CY leads to a discrete approximation Ah to 
A. The finite element approximation to ([21) is then given by 

(3) A'^^TTh + Ah)~^^Thd^i, 

^ Jo 

where tt/j is the L^(D)-projection onto Yh- In [Ilj . we proved the convergence in 
of an equivalent version of this method when A = S' is real, symmetric and 
positive definite. We also showed the exponential convergence of a sine quadrature 
approximation. 

The current paper extends the approach of m to the case when A is a regularly 
accretive operator. The proof provided in m is based on the fact that, in the 
Hermitian case, the domain of S^, for 7 € R, is naturally characterized in terms 
the decay of the coefficients in expansions involving the eigenvectors of S. Assuming 
elliptic regularity, it is then possible show that D(S®^^) for 0 < s < I -I- a coincides 
with standard Hilbert spaces. Here a is the regularity parameter (see below). Thus, 
norms of the operator {fil + S)~^ acting between the standard Sobolev spaces 
can be bounded using their series expansions, the norms in and Young’s 
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inequalities. In contrast, the spaces D{A^) cannot be characterized in such a simple 
way when A is not Hermitian. 

The main result of this paper is the following error estimate (Theorem 16.21 and 
Remark IQ) : for any 0 < r < 1 there exists <5 > 0 and a constant C independent 
of h such that for / G D(A^) 

( 4 ) UA-P - A-^7Th)f\\Hr(,n) < 

with C being replaced by C'log(h“^) for certain combinations of r, a, S and /3. 
Here a > 0 is the so-called elliptic regularity pick-up which is the regularity above 
expected for A~^ f for appropriate /, see Assumption [TJ Even in the case 
of Hermitian A, the above result extends those in m to r > 0. 

This paper shows that the general approach for proving (jd]) in m can be ex¬ 
tended to the case of regularly accretive A with additional technical machinery. 
Some of the most challenging issues involve the relationship between for 

s G [0,1 -I- a] and fractional Sobolev spaces. The case of s G [0,1) is contained in 
the acclaimed paper by Kato m showing that for regularly accretive operators, 
£)(^s/ 2) coincides with the interpolation space between L^(n) and V defined using 
the real method. The case of s = 1 is the celebrated Kato Square Root Problem. 
This is a deep result which has been intensively studied (see, [niiisniET] and 
the references in [39]). The results in those papers give conditions when one can 
conclude that D{A^/^) = ZI((A*)^/^) = V, where A* stands for the adjoint of A. 
Motivated by the approach of Agranovich and Selitskii [1] for proving the Kato 
Square Root Problem, we show in this paper that under elliptic regularity assump¬ 
tions, iJ®(r2) n V C ZI(A®/^) and i7®(H) fl V C Z3((A*)®/^) for s G [0,1 -I- a] with 
equality when additional injectivity assumptions on A~^ and (A*)“^ hold. With 
this information, the norms of {fil -I- acting between Sobolev spaces can be 

bounded in terms of the L^(H) operator norm of A*(/r/ -|- A)~^ with t G [0,1]. 
This, in turn, (see Lemma l6.3|l can be bounded by interpolation using the fact that 
D{A^) coincides with the interpolation scale between L^(n) and D{A) (using the 
complex method). 

Similar results for 0 < s < 1 are also required for the finite element approxi¬ 
mation Ah- In a discrete setting, the question is to guarantee the existence of a 
constant C independent of h such that for all Vh G Yh 

This is provided by Lemma 15.21 for 0 < s < 1 and a solution to the discrete Kato 
problem (s = 1) is given in Theorem 16.51 

We also study a sine quadrature approximation to AJ^^f for / G V^. A change 
of integration variable shows that 

(5) A-^ = r e(i-«^(e^/ + A,)-^ dy. 

^ J — OO 

Motivated by [36], the sine quadrature approximation to AJ^^ is given by 

(6) Q-/(A,) := •£ + A,r'. 

TT 

i^-N 

Here k := I/Vn is the quadrature step size, A^ is a positive integer and ye := £k for 
£ = —N, The standard tools related to the sine quadrature together with the 
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characterization of for s S [0,1/2] mentioned above yield the quadrature 

error estimate iRemark 17.21) 

where Cq is a constant independent of k and h. 

The outline of this paper is as follows. In Section [21 we introduce the notations 
and properties related to operator calculus with non-Hermitian operators. Sec¬ 
tion [3] is devoted to the study of the Hermitian part of A and the related dotted 
spaces. The latter is instrumental for the characterization of D{A'^) discussed in 
Section [d] The finite element approximations are then introduced in Section |5l 
which also contains the proof of the error estimate (|4]). This coupled with the ex¬ 
ponentially convergent sine quadrature studied in Section [721 yields the final error 
estimate for the fully discrete and implementable approximation. We end this work 
with Section [5] providing a numerical illustration of the approximation of fractional 
convection-diffusion problems. 

2. Fractional Powers of non-Hermitian Operators 

We recall that H is a bounded domain in with a Lipschitz continuous boundary 
r which is the disjoint union of an open set T^v and its complement To = T \ T^v- 
Let Lp'iTl) be the space of complex valued functions on H with square integrable 
absolute value and denote by ||.|| and (.,.) the corresponding norm and Hermitan 
inner product. Let be the Sobolev space of complex valued functions on Tl 

and set 

V := {u G u = 0 on r£)}, 

the restriction of H^{Tl) to functions with vanishing traces on Tu. We implicitly 
assume that L^i is such that the trace operator from H^{Tl) is bounded into L^(To), 
e.g., r D does not contain any isolated sets of zero d—2 dimensional measure. We 
denote ||.||i and ||.||v to be the norms on H^{Tl) and V defined respectively by 



For convenience, we avoid the situation where the variational space requires L^(H)- 
orthogonalization, i.e., the Neumann problem without a zeroth order term. 

For a bounded operator G : X{Tl) — Y{Tl) between two Banach spaces {X (H) ,IMU) 
and (F(H), ||.||v) we write 

||G||x^y := sup ||Gu||v 

u^X,\\u\\x = l 

and, in short, ||G|| := ||G||z, 2 ^i 2 . 

Throughout this paper we use the notation A < B to denote A < CB with a 
constant G independent of A, B and the discretization mesh parameter h (defined 
later). When appropriate, we shall be more explicit on the dependence of G. 

Consider the sesquilinear form o for u, V in V. We assume that A{-, •) is strongly 
elliptic and bounded. The assumption of strong ellipticity is the existence of a 
positive constant cq satisfying 

(7) lHe(H(i;,'u)) > co||u||v. 


for all u G V. 
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The boundedness of on V implies the existence of a positive constant a 

satisfying 


(8) |yl(u,w)| < ci||u||v||u||v, forallu,ueV. 

The conditions 0 and ([U imply that the sesquilinear form A(-, •) is regular on V 
(see, Section 2 of IS]). 

Following m, we define the Hermitian forms 


^ , A(u,v) + A(v,u) , ^ , 

D\eA{u,v) := - - - and JmA{u,v) := 

Note that o implies that D\eA{u,u) is equivalent to ||u|ly, 
equivalent to 


A{u, v) — A{v, u) 

2i ■ 

for all u S V and ([5]) is 


(9) \3mA{u,u)\ < T]9\tA{u,u), for all u G V, 

for some rj > 0. The smallest constant j] above is called the index of A{-^ •). 

We now define operators associated with regular sesquilinear forms. Let A be 
the unique closed m(aximal)-accretive operator of Theorem 2.1 of |31] . which is 
defined as follows. We set T : L'^iVL) —>• V by 


A{Tu, (j)) = (u, (j)), for all (/) G V 


(which is uniquely defined by the Lax-Milgram Theorem) and set 

D{A) ■- Range(f) c V. 

As T is one to one, we define Aw := T~^w, for w G D{A). The operator A 
associated with a regular sesquilinear form is said to be regularly accretive. 

We denote by V* denotes the set of bounded antilinear functionals on V. It 
will be useful to consider also a related bounded operator : V* —?> V defined for 
F G V* by the unique solution (Lax-Milgram again) to 

A{TaF,(l)) = {F,4>), forall<(.GV. 

Here (•,•) denotes the antilinear functional/function pairing. The operator Ta is a 
bijection of V* onto V and we denote its inverse by A^. ^From the definition of Tq, 
we readily deduce that for u G V, AaU satisfies 


( 10 ) 


{AaU, v) = A{u, v), for all u G V. 


It is also clear from the definition of A that u G D{A) if and only if AaU extends 
to a bounded antilinear functional on L^(H) and, then. 


{Au,v) = {AaU,v) = A{u,v), for all u G D{A),v G V. 


The above constructions can be repeated for adjoints defining V* the set of 
bounded linear functionals on V, T* : L^(n) —)> V, D{A*) := Range(T*) C V, 
A* := (r*)-i : ^(A*) ^ T* : V; ^ V and A* := {T*)-^. In this case, for 

u G V, Aiv satisfies 


(11) {u,Aiv)=A{u,v), for all u G V. 

with (.,.) also denoting the function/linear functional pairing. We also have that 
V G D{A*) if and only if Aiv extends to a bounded linear functional on L^(fl) and, 
then. 


{u,A*v) = {u,A*v) = A{u,v), for all u G V, u G D{A*). 
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Of course, these definitions imply that for u S D{A) and v S D{A*), {Au,v) = 
(u, A*v). 

By construction, the operators A and A* defined from regular sesquilinear forms 
are regularly accretive (cf. [31]). They satisfy the following theorem (see, also [23]): 

Theorem 2.1 (Theorem 2.2 of [31]). Let A he the unique regularly accretive opera¬ 
tor defined from a regular sesquilinear form A(-, •) with index g. Set uj := arctan(77). 
Then the numerical range and the spectrum of A are subsets of the sector S^j ■= 
{z G C : I argz| < cu}. Further, the resolvent set p{A) of A contains 50 :=C\Scj 
and on this set the resolvent Rz{A) := {A — zl)~^ satisfies 

IIP MJII <-/ w < arg(z) < f + w, 

llltzfAjll S I |z|-i arg(z)>f+w. 

The result also holds for A replaced by A*. 

Remark 2.1. It easily follows from 0 that for fHe(z) < co/2, 

y||u||? < \Aiu,u) - z{u,u)\ 

which implies that 

||i?.(A)/||i< -ll/ll. 

Co 

It follows from Theorem [2T] and Remark [2T] that the Bochner integral appearing 
in for fi G (0,1), is well defined and gives a bounded operator A ^ on Lflfil). 
Fractional powers for positive indices can be defined from those with negative in¬ 
dices. For fi G (0,1), 

D{A^) = {u&L^{n) : A^-^u&D{A)} 
and A^u := A{A^~^)u for u G D{A^). 

An alternative but equivalent definition of fractional powers of positive operators 
(for fi G (—1,1)) is given in, e.g., [35]. We shall recall some additional properties 
provided there (for fi G (0,1)). Theorem 4.6 of [35] implies that D{A) C D{A^) 
and for v G D{A), A^v = and Av = Also, for any 

/3 > 0, D{A^) = {A~^v : v G L^(fl)} and A^v = {A~^)~^v for v G D{A^). 

Set w := A^v. The last statement in the previous paragraph implies that w = 
{A~^)~^v. Now as {A~^)w = v, 

IHI = p-''^c||<||A-/^||||u;|| = ||A-/5||||A^||. 

This implies that we can take 

(12) |k||D(Aq := 

as our norm on D{A^) for s G [0,1). 

Using the above and techniques from functional calculus [26] . we can conclude 
similar facts concerning products of fractional powers and the resolvent, namely, 


(13) 

and 

A-f^ Rz{A)u = Rz{A)A-^u, 

for all u G Lf{FL),j3 > 0, 

(14) 

A>^Rz{A)u = Rz{A)A^u, 

for all u G D{A^),P G [0,1], 


We shall also connect fractional powers of operators with their adjoints in the 
L^(n)-inner product. We have already noted that for u G D{A) and v G D{A*), 
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{Au,v) = {u,A*v). This holds for fractional powers as well, i.e., {A^u,v) = 
{u, {A*)^v) provided that u S D{A) and v S D{A*). 


3. The Hermitian Operator and the Dotted Spaces. 

For notational simplicity, we set S{u, v) := v). As already noted, S{u, 

provides an equivalent norm on V and we redefine ||u||v := S{u,uy^‘^. As 
is regular (i.e. satisfies m and with A(.,.) replaced by S {.,.)), there is an as¬ 
sociated (m-accretive) unbounded operator S. The latter is defined similarly as A 
from A(.,.) in Section [5] (see also [3T]). This is, upon first defining Ts : T^(0) —>• V 
by Tsf '■= w where u> £ V is the unique solution of 

S{w, (j)) = (/, <^) for all (() G V 

and then setting D{S) := Range(rs), S := Tg^. 

In addition, as is symmetric and coercive, S is self-adjoint and satishes 

S{u,v) = S^^^v). 

We consider the Hilbert scale of spaces defined by := D(S'®) for s > 0. The 
above discussion implies 


(15) 


=■ 


and 




Moreover, the operator Ts is a compact Hermitian operator on L^(n) and so there 
is a countable L^(H)-orthonormal basis {ipi, z = 1,..., oo} of eigenfunctions for Tg. 
The corresponding eigenvalues {/Zi} can be ordered so that they are non-increasing 
with limit 0 and we set Ai = This leads to a realization of in terms of 

eigenfunction expansions, namely, for s G (0,1) 

f OO OO N 

(16) 11(5'®) := iil® = i w = G L^(H) : ^ |(w,'(/'i)pA® < oo 1. 

^ i=i j=i ^ 

The spaces FI® are Hilbert spaces with inner product 

OO 

i=i 

Moreover, they are a Hilbert scale of spaces and are also connected by the real 
interpolation method. 

As already mentioned = V so that the set of antilinear functionals on V, 
denoted V*, can be characterized by 


ra=Ha^ 


E 

i=i 


^A, ^ < 00 




where J2T=i CjV'j, j ■= i Cjdj. 

In addition, the set of antilinear functionals on L^(H), denoted by L^(H)*, is 


given by 


f j OO y OO 

^ ' 7=1 ' 7=1 


< (X) 
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Hence, the intermediate spaces are defined by 


:= 



OO 

El 

i=i 


12 \ —s 


A"'* < OO 


and are Hilbert spaces with the obvious inner product. These also are a Hilbert 
scale of interpolation spaces for s S [—1,0]. In addition, these spaces are dual to 
H^, i.e., if s G [0,1] and 


then 

(17) 


(w,-) 



G H- 


W\\^-s = 


OO 


1=1 


1/2 


= sup 


{w,0) 


and if 0 G i7®, 
(18) 




sup 


{w,0) 


Considering linear functionals instead of antilinear functionals and replacing V* 
and L^(H)* by spaces of linear functionals, V* and L^(fl)* gives rise to the analo¬ 
gous Hilbert scale with = V* and = L^(n)* as endpoints with equalities 
similar to (HU and (fTSl) holding for these as well. 

As we shall see in SectionSJ D{A^/^) relates either to iJ® or 77®(n)nV depending 
on whether s > 0 is smaller or greater that 1. In order to unify the presentation, 
we introduce the following spaces equipped with their natural norms: 


:= 


77®(fl)nv 


for s G [0,1], 
for s > 1. 


4. Characterization of D{A^ 

In this section, we first observe that the dotted spaces for s G [0,1) coincide with 
the domains of fractional powers of A and A* (c.f., m)- In addition, we note that 
the dotted spaces and can be identified with the dual space of iJ®, for 
s G [0, 1]. The case of i7® := i7®, s G (0, 1) is addressed in the following theorem 
which is an immediate consequence of Theorem 3.1 of Ell- 

Theorem 4.1 (Characterization of D{A^) for 0 < s < 1). Assume that ([7]) and 
(ID) hold. Then for s G [0,1), 

ZI(A®/2) = Z1((H*)®/2) = D(5'®/2) = H®, 

with equivalent norms. 

The identification of the negative dotted spaces with with the duals is given in 
the following remark. 

Remark 4.1 (Characterization of Negative Spaces). We identify f G Lf{VL) with 
the functional Ff G L^(H)* defined by 

{Fi,9) = {f,0), for all0&L\n). 
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It follows from Theorem \4-l\ and (1171) that the norms and ||A '’/|| are 

equivalent (for s € [0,1/2 )). Indeed, 


\\Pi\\H7^‘= sup 




{fL4>) 


sup 


if, </) 


= sup 

&GL^{n) 


if, iA*)-^e) 


p-vii- 


Here ~ denotes comparability with constants independent of f. For simplicity, we 
shall write ||/||^- 2 s instead of ||F/||^- 2 s. We can identify with L^{Tl)( in 

an analogous way and similar norm equivalences hold. 


Elliptic regularity is required to obtain convergence rates for finite element ap¬ 
proximation. Such results for boundary value problems have been studied by many 
authors 0 in uni 123131131 isi sni m] ■ The follow assumption illustrates the type 
of elliptic regularity results available. 

Assumption 1 (Elliptic Regularity). We shall assume elliptic regularity for the 
form A(-,-) with indices a € (0,1]. Specifically, we assume that for s € (0, a], Ta 
is a bounded map of into and T* is a bounded map of Hi into 

H^+F 


The above assumptions imply the following theorem. 

Theorem 4.2 (Property of D{Ai) for s > 1). Assume that (l7|), ([8]) and the elliptic 
regularity assumptions (Assumption]^ hold. Then for s G (1,1-1- a], 

D{A‘‘/^) c H^ and D{{A*y/^) c Hf 
with continuous embeddings. 

Remark 4.2 (Kato Square Root Problem). The case of s = 1, i.e., D{A^^^) C 
V =: H^ with continuous imbedding is contained in the Kato Square Root Theorem. 
This is a deep theorem which has been intensively studied, see [I E El 132] and 
the references in |39j . The Kato Square Root Theorem holds for our problem under 
fairly weak regularity assumptions on the coefficients defining our bilinear form [1] . 
In fact, it requires the existence of e > 0 such that Aa and A/ are bounded operators 
from to , respectively, for jyj < e. 

Proof of Theorem \4.2\ We consider the case of A as the case of A* is similar. Sup¬ 
pose that u is in D{A) and v is in D{A*). Then, for t G [0, a], 

A(u,v) = (Au,v) = 

Thus, Theorem Q] gives 

\F{v)\ < ||A(l+*)/2u||||(^*)(l-t)/2^|| 

This implies that F G . The elliptic regularity Assumption [T] implies that 

u = TfF is in and satisfies 

As D[A) is dense in D{A^^+*F^), £)(A(i+*)/2) c H^+* follows. □ 
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5. Finite element approximation 


In this section, we define finite element approximations to the operator A~^ for 
/3 € (0,1). For simplicity, we assume that the domain is polyhedral so that 
it can be partitioned into a conforming subdivision made of simplices. Further, 
we assume that we are given a finite dimensional subspace C V consisting of 
continuous complex valued functions, vanishing on F^i, which are piecewise linear 
with respect to a conforming subdivision of simplicies of maximal size diameter 
h < \. Notice that when the form A(v,w) is real for real v,w, so is the finite 
element space (see Remark 15.11) . We also need to assume that the triangulation 
matches the partitioning F = F^i U F^v. This means that any mesh simplex of 
dimension less than d which lies on F is contained in either Fjv or Fd. Given a 
universal constant p > 0, we restrict further our considerations to quasi-uniform 
partitions T, i.e. satisfying 


(19) 


maxTgT- diam(T) 
minj'gT- diam(T) 


< P- 


Let TTh denote the L^(n)-orthogonal projector onto Yh- Given a sequence of 
conforming subdivisions {T}h satisfying (fT^ . there holds 


( 20 ) 


hhv\\i < Cllull 


for all u G V, 


where the constant C is independent of h] see [M]. Obviously, nh is a bounded 
operator on L^(fl) and, by interpolation, is a bounded operator on iJ® for s G [0,1] 
with bounds independent of h. 

We shall need the following lemma providing approximation properties for tt/i. 
These results are a consequence of the Scott-Zhang approximation operator |43j 
and operator interpolation. As the specific form of these results are needed for the 
analysis in the remainder of this paper, we include a proof for completeness. 

Lemma 5.1. Let s be in [0,1] and a > 0 be such that s a < 2. Then there is a 
constant C = C{s,a) not depending on h and satisfying 


\\{I-TTh)u\\fj, < Ch‘^\\u\\ 




for all u G iL®+'". 


Proof. Let Wh denote the Scott-Zhang approximation operator [43] mapping onto 
the set of piecewise linear polynomials with respect to the above triangulation 
(without any imposed boundary conditions). This operator satisfies, for £ = 0,1 
and fc = 1,2, 


( 21 ) 


\\{I - Trh)u\\He A h^ \\u\\h''‘, for all u G iF (O). 


In addition, iThU G Yh for u G V. 

We first verify the lemma when s -|- cr < 1. We clearly have 

(22) ||(J - 7r/i)u|| < 11(7 - 7f/i)M|| ^ /i||m||v, for all m G 
It immediately follows from (1^01) that 

(23) II(J — 7r/i)M||v ^ ||m||v, for all u G V. 
Interpolating this and (|^ gives 


(24) 


\\{I -TTh)u\\pj, A 1|u||v, 


for all u G V. 
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for s S [0,1]. As TTh is stable on and V, interpolation implies that it is stable 

on 77® and hence 

(25) IK-f-^ llulljj., for all M e ff®, 

for s G [0,1]. Interpolating (l24l) and (I25|) and applying the reiteration theorem 
gives for cr > 0 and s + cr < 1, 

(26) ||(/-^i.)u|b. for all u e 77*+^ 

We next consider the case when s + a G (1,2]. Taking 7 = 1 in (|2T]) and 
interpolating between the fc = 1 and 7 = 2 gives for ( G [1,2], 

11(7 - nh)u\\Hi ^ h‘^~^\\u\\HC, for all u G 77^(71). 

Thus for u G ff®+'^, by ([261), 

\\{I - Trh)u\\fj, di /i^“®||(7- 7r/i)u||v ^ /i^“®(||(7 - ^h)u||v + U^h - 7r/i)u||v) 

^ +/i-®||(T?,-7r?,)u|| ^ 7'"|[u||jj,+,, 

where the last inequality followed from (EU and obvious manipulations. □ 


We define Ah :Yh ^ V/, by 

(27) {AhVh,iph) = A{vh,iph), ioT all (fh GYh- 

The operator Ah is the discrete analogue of A and we analogously define A)), the 
discrete analogue of A*. The fractional powers AJ^^ for /3 > 0 are again given by 
@ but with A replaced by Ah, i.e., for /3 S (0,1), '-Yh ^Yh is given by 

(28) A-^ := r + Ah)-^ dfi. 

TI" Jo 

The goal of this paper is to analyze the error between A~^ f and A^^nhf. 

Remark 5.1 (Real Valued Bilinear Forms and Finite Element Spaces). When 
A(v,w) is real for real v,w, the above operators restricted to real valued functions 
are real valued and hence we may use Sobolev spaces and approximation spaces Yh 
of real valued functions. 


Similarly, let Sh '-Yh ^ Yh be defined by 


{ShVh,Wh) = S{vh,Wh), for all Vh,Wh G Yh- 


Theorem 3.1 of [3T] applied to the discrete operators Ah and Sh shows that for 

se [0,1/2), 


1 — tan 


TTS 


||^^/.|| < ||A^u,|[ < 


1 + (— tan7rs)^^^(77 + 


II 


for all Vh G Yh- Here rj is the index of A(-, •) (see, dH])). This also holds for A)). 

The bound (1^ implies (see, e.g., [7]) that there are positive constants c and C, 
not depending on h such that for s G [0,1], 

(29) c\\Sl^^Vh\\ < ||p/.||^. < C\\Sl^\h\\, for all Vh G Yh- 

Combining the preceding two sets of inequalities proves the following lemma. 
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Lemma 5.2 (Discrete Characterization of for s G [0,1)). There exist positive 
constants c, C independent of h such that for all Vh G Yh and s G [0,1), 

cWvhWns < \\A"f/'^Vh\\ < C||p;,||^.. 

This result holds with Ah replaced by A'^. 

Remark 5.2 (Quasi-uniformity Assumption). The above lemma still holds as long 
as (1201) holds. This is for instance the case for certain mesh refinement strategies 

BM- 


6. Error Estimates 

In this section, we study numerical approximation to the operators A~^ for 
/3 G (0,1). Specifically, this involves bounding the errors {A~^ — Af^^TTh)v for v 
having appropriate smoothness. 

We shall use our finite element spaces to approximate A~^ := T^. Specifically, 
Th,a : V* -)■ V/j is defined for F e V* by 

A{Th,aF, (fh) = {F, fih), for all 4>h G Yh- 

We define T^ i corresponding to T* := (A*)“^ analogously. 

The following lemma provides approximation error bounds in terms of the norms 
needed for our subsequent analysis. Although the techniques in the proof (Galerkin 
orthogonality and Nitsche finite element duality) are completely classical, the re¬ 
sults are not quotable (as far as we know). We include a proof for completeness. 

Lemma 6.1 (Finite Element Error). Assume that (jT]), (0 and the elliptic regularity 
AssumptionUi hold. Let s G [0, |] and set a* := ^{a + min(l — 2s, a)). There is a 
positive constant c not depending on h satisfying 

(30) \\Ta - Ta,h\\HS-^^H^s < ch^^\ 

The above immediately implies 

\\Ta-Ta,h\\<ch^'^. 


Proof. The proof of this lemma is classical and we only include details for com¬ 
pleteness. We distinguish two cases. 

□ When 2s < 1 — a then we can fully take advantage of the elliptic regularity 
assumption. For F G we set e = (Tq — Ta^h)F. By (fTSl) and the elliptic 

regularity Assumption [1] 


< ||e||^i-„ A sup 
G&Hi 


(e,G) 


- I|G|| 


A sup 


A{e,TfG) 




= sup 


A(e, w) 


A{e,w- Wh) 


= sup 


\w\\hi+^ 


We used Galerkin orthogonality for the last equality (which holds for any Wh G Yh). 
Using the approximation property 

inf ||iT — ru/illi ^/i“||u;||^i+c, for all u; G 
wheYh 

0 and the above inequalities give 
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[J] This duality argument yields a reduced order of convergence when 2s > 1 — a. 
Indeed, proceeding similarly 

„ „ ^ (e,G) ^ Aie,T;G) 

lle||^ 2 . ^ sup —- ^ sup ——- 

||G||^-2, ceH-^‘ Pi ^\\h2-2s 

A(e,w) A(e,w — Wh) 

= sup -= sup —-, 

u;gLr2-2s ||'Ui|| ff2-2s j„g//2-2s 11 "Ui 11 ff 2-2s 

SO that together with the approximation property 

inf ||w — w/i||i ^/i^“^®||u;||// 2 - 2 s, for all ic S 

uiheVh 


we deduce that 


|e||^2. ^ 


Gathering the two cases 2s > 1 + a and 2s < 1 + a, we get 

Whence, together with the estimate 

WehAh^^WTaFWHi+.Ah^FW^.-^ 

guaranteed by Cea’s Lemma and elliptic regularity Assumption [TJ we obtain (I5D1) 
as desired. □ 


We can now state and prove our main convergence results. It requires data in 
the abstract space D[A^) for some 5 > 0. A characterization of D{A^) is provided 
in Theorem 16.41 below. 

Theorem 6.2 (Convergence). Suppose that ([7|) and ([5]) as well as the elliptic 
regularity Assumption\^hold. Given s S [0, ^), set a* := i(Q;+ min(l — 2s, a)) and 
7 := max(s + a* — /3,0) and let ^ > 7 . Assume finally, that s + a* 1/2 when 
(5 = 7 and s + a* > /3. Then there exists a constant G independent of h and S such 
that 

UA-P-Af7Th)f\\H2s<Cs,hh^^^\\A^f\\, for all f € DiA^), 

where 

{ G ln(2 / h) : when ^ = 7 and s + a*>/3, s + a^, ^ — 

G : when 5 > 7 and s + a^, > fi, 

G : when 6 = 0 and /3 > s + Of*. 

Remark 6.1 (Critical Case s + a* = ^). The condition s + a* ^ in the above 

theorem can be removed provided that the Kato Square Root Theorem holds as well 
(see, Remark \J^ . 

Remark 6.2 (Critical Case 2s = 1). The above results also hold when 2s = 1 
provided that the continuous and discrete Kato Square Root Theorem hold. As 
already mentioned above, the former relies on the additional assumption requiring 
the existence of e > 0 such that Aa and A( are bounded from F[^~^'^ to H2~^ 
and , respectively, for I 7 I < e. For the discrete Kato Theorem, we will need 
to assume similar conditions for operators based on the S form, see Theorem | 6 ‘./| 
below. 


The above theorem depends on an auxiliary lemma. 
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Lemma 6.3. For s G [0,1], there is a constant C not depending on h such that for 
any fi G ( 0 , oo) 

\\A‘^{g, + for all v € L'^{Vl) 


and 


\\Al{p + Ah) Vh\\<Cfj.'‘ ||i;/i||, for allvhGYh- 


Proof. The claim relies on interpolation estimates. As the same argument is used 
for both estimates, we only prove the first. 

Theorem 4.29 of [3S] implies A** is a bounded operator satisfying 

(32) ||A*‘|| < forallteR. 

This, in turn, implies that (e.g., Lemma 4.31 of [35]) for s G [0, ^], 

(33) [L\n),D{A^/^)]2s=DiA^). 

Here [X, Y] 2 s denotes the intermediate space between X and Y obtained by the 
complex interpolation method. Thus, for w G D{A), Corollary 2.8 of [S^ gives 

IIA^^II A |k||[L^(0),D(A)]. < ||An;|n|R;|^-^ 

Now ii w = {pL + A)~^V with v G then 

/i||w|p + lHeA(w, w) = w) 

and hence 0 immediately implies ||(/i +A) ^v\\ < p ^||n||. In addition, 

||Aw|| = ||i; - p{p + A)-in|| < 2||i;||. 

The lemma follows combining the above estimates. □ 

Proof of Theorem 1 6. ‘A Without loss of generality, we may assume that (5 < 1 + a, 
since we shall get 2 q!* order convergence as soon as i5 > 2a, — 2/3 and we always 
have 2 a, — 2,5 < 1 + a,. 
nn We first show that 

(34) \\{I-^n)A-Pf\\p,.Fh‘^<^‘\\A^f\\. 

Theorem 4.6 of |3Sj implies that A“^/ is in D{A°^‘’) when / is in D{A°‘*~^) and we 
now discuss separately the cases s+a, G (0, 5 ), s+a, = ^ and s+a, G (^,^( 1 + 0 )]. 
When s + a, G (0, 4), we apply Theorem 14.II and obtain 


||(/ - 7r^,)A-Pf\\p,, + IIA^+^-'^/ll 

+ h 2 “*||A^/||, 

recalling that 5 > 7 >s + a, — /3. For s + a, G (^, 4(1 +a)], we apply Theorem l4.2l 
to conclude that A~^f is in iy 2 s+ 2 a, ^ again, 

+ h 2 “*||AV|l. 

Finally, we consider the case s + a, = which entails a < 1 — s so that a, = a 
and 2s + 2a = 1. We choose 0 < e < a (further restricted below) so that as above 

11(7 - nh)A-Pf\\p,. + h^^+^\\A-^ + h2“+^||A^-/5/||. 
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In addition the assumption ^ > 7 := max(i — /3,0) yields 5 “ /? + f < ^ upon 
choosing a sufficiently small e. Hence, we deduce 

This, (1551) and (I5S1) yield (IMl) . 

[ 2 ] By the triangle inequality, it suffices now to bound 


(37) 




< 


sin(/37r) 




\TTh{ll + A) ^-{^l + Ah) d^l. 


Assuming without loss of generality that /i < 1, we shall break the above integral 
into integrals on three subintervals, namely, ( 0 , 1 ), and 00 ). 

0 We start with 00 ). Recalling the definition of the operator norm (fT^ 

as well as the characterizations of the dotted space provided by Theorem 14.11 , we 
get 

where we used in addition the stability of TTh in D{A^) and the boundedness of A“’’ 
from L^(H) to for r > 0 (see discussion below Remark l2.1|) . Hence, applying 

Lemma 15.21 yields 

Similarly, but using the discrete characterization provided by Lemma l5.21 we obtain 
||(ai +^ ^max(.-yo)-l^ 

Thus, invoking Lemma [6.31 we deduce that 

pOO 

/ ^\\'^hi^J- + A) ^ — {fj, + Ah) 




Jh 


!c( —a*, —/?) — 1 


h - 2 c ^*/0 


dfi :< 


because S > s + a^: — /3. 

0 For the remaining two cases, we use the identity 

'^/t(/^ + A) ^ — (fj, + Ah) ^TTh = {fJi + Ah) ^AhT^hiTa ~ Th,a)A{fl + A) 


The latter, follows from the identification of Remark |4T] and that the observation 
that for u £ D{A) and u G V, 

(TaAujv) = A{TaAu,T*v) = {Au,T*v) = A{u,T*v) = {u,v), 
i.e., TaAu = u. Also, it is easy to see that AhnhTh,a = T^h- Thus, for u G D{A), 
AhTTh{Ta - Th,a)Au = + Ah)TThU - TT/j (/i + A)u, 


which leads to the desired identity. 
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[H For fi G (1, 2 a.//3^^ write 

I M ^ ^h) ^T^hiTa — Tii^a)A{ll + A) ^\\dfl 

M ^\\Ah{lJ^ + Ah) ^7r;i||^i_c._j^^2s ||(T'a ~ A ,j)||^-l+<»,_^^i_„^ 

||24(/r + A) ||£)(2i'5)^ij~'-+“* 

Now, the definition (1121) of ||■||D(A'5) together with the characterization of the 
negative spaces provided in Remark 14.11 imply that 

In addition, using Lemma 15.21 we obtain 

U A !■ , A \-l II _ o \\Ah{^J' + Ah) ^TThfll^^s 

\\Ah{tJ' + Ah,) TT^ II jjl-o:, — sup 

||^/t(/^ F 'n'hf \\]J 2 a 


:< sup 
fem- 




IIA " '^\t^ + Ah) I5II 


d: sup 
ge^Vh 
-1+' 


-+S 


The above three estimates together with Lemma I5TT] yield 

M ^'^Ah{^-\-Ah) ^"KhiTa ~ Th,^a)A{^-\-A) jj(^jys^^fj2a dfj, 




-l+s+a.-P-S j 

M ^2a. 


^ 2 a, \TiY[2/h) : if (5 = s + q;. — /3, 

: otherwise. 


0 Finally for /r £ (0,1), we write 

f /r ^\\Ah{f^-\-Ah) ^TTh{Ta ~ Th,a)A{fi-\-A) dfi 

0 

— / M ^||^/t(li + 4l?i) 7r/i||^2s_j,^2s ||7a — T;i^a||j^2_j,^2s ||^(/r + A) ||l 2_,,^2 d/r 

Jo 

d / ij,~^ dfi d hd°‘. 

Jo 

0 The proof of the theorem is complete upon collecting the estimates obtained in 
Steps 3,5 and 6. □ 


Theorem Q] and the Kato Square Root Theorem characterize D{A^) for s £ 
[0,1/2]. The characterization can be extended to s £ (1/2, (1 + a)/2] when Aa 
maps into This is of particular importance to characterize the regularity 

assumption / £ D{A^) in Theorem l6.2l 

Theorem 6.4 (Characterization of D{AA+^'>/'^) for s £ (0,a]). Suppose that © 
and ([8]) hold. Assume furthermore that for s £ (0,a], 
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(38) Ta is an isomorphism from ^^'^into 
Then, 

with equivalent norms. 

Proof. By Theorem 14.21 we need only prove that C We first 

observe that D{A) fl is dense in Indeed, if w is in then (l3^ 

implies that AaW is in . As is dense in there is a sequence Fn S 

converging to AaW in Setting Un '.= {Aa)~^Fn, elliptic regularity implies 

that Un converges to w in Clearly u„ is in D{A), i.e., D{A) fl is dense 

in Rl + S 

as claimed. 

Suppose that u G D{A) n We first show that 

(39) 

For V G D{A*) and 5 := (1 - s)/2 G [0,1/2), 

(yl(i+-)/2y^^) = {A-^Au,v) = (Au, iA*)-^v). 

Since u G D{A) and {A*)-^v G D{A*) c V, 

|(Azi,(A*)-M| = \{AaU,{A*)-'v)\ < |lA„n||^-.+.||(A*)-^||^,_. 

where we also used CZD and ((38l) . Now, Theorem 14.11 ensures that 
||(A*)-^i;||^,_. ^||(A*)^(A*)-^u|| = ||t||. 

Combining the above inequalities shows that (IM)) holds for u G D{A) fl 77^+®. The 
inclusion C Z1(A(^+®)/^) and (l3^ for u G hold by density. □ 

Let Ts,a be defined similarly to Ta but using the form £'(•,•). The final result 
in this section shows that under suitable assumptions. Lemma |5.2I holds for s = 1. 
This is a discrete Kato Square Root Theorem. Its proof was motivated by the proof 
of the Kato Square Root Theorem given in [1] . 

Theorem 6.5 (Discrete Kato Square Root Theorem). Suppose that 0 and ([8|) 
hold. Assume furthermore that for some s G (0,1/2), Ta, T* and Ts,a are isomor¬ 
phisms from into iL^+®. Then, there are positive constants c,C independent 

of h such that 

c||nft||v < W^H'^VhW < CIIu/jIIv, for all vu G Vh- 
The analogous inequalities hold with A’^ replacing Ah above. 

Proof. In this proof C denotes a generic constant independent of h. The assumption 
on Ts^a implies that for t G (0, s], with equivalent norms [llj . 

We first observe that it suffices to show that 

(40) \\A]!\h\\ < C||iq.||v, for all Uh G V„ 

along with the analogous inequality involving Indeed, if (l40)) holds then by (O, 
colk.lll < \{AhUh,Uh)\ < \\A^'\hmAlY/^uh\\ 

< C'||(A)))^/^u/j||||m/i||v 
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and hence 

coWuHh < C\\iAl)^/W\\. 

1 /2 

The proof of the lower bound involving is similar. 
Applying Lemma gives, for Uh S Vh, 

(41) 


W^h “/ill = sup 

^ C||< C'||m;i||^i+s 


se^h 


||(A*)(i-)/2d|| 


{Auh,0) 

- ^ - 

eeNh 


where we used the assumption on Ta = (Aa)“^. 

Let TTh denote the S'-elliptic projection onto Yt, i-e., for v € Y, Wh = ^hV € Yh 
solves 

S{wh,0h) = S{v,0h), for all 0,1 e V/j. 

It is a consequence of the isomorphism assumption on Ts^a that tt/j is a uniformly 
(independent of h) bounded operator on (see, e.g., [13] )• Thus, recalling the 

eigenvalue decomposition (HSj), it holds for Uh & Y^, 

Siuhjfj)) S{uh,'^h(l))\\^h(t>\\H^+. 

WUhllHl-s = sup -y—- = sup ——-- 

, ^ S{uh, iph) 

- ^ ^'^p Tin—■ 


Similarly, for all m/i G V/i, 


|u?i||^i+. = sup 


S{uh, 4>h) 


fjl-B 


Thus, the characterization (l29l) gives 

\\ndH^.B<C sup = 

4>h€V^ IIS-^ 

Combining this with (IdTI) gives 

for all uh S V„. 

Interpolating this result with the trivial inequality 

||A°u/i|| < ||S'°Uft||, for all e V/i, 

gives 

WAl^^UhW < C\\Sl^\h\\ = C||u,,|lv, for all S V,, 
and verifies (ini)- The proof for the analogous inequality involving A'^ is similar. □ 


7. Exponentially Convergent sing Quadrature. 

Theorem l6.2l and Remark ldTTl provide estimates for the errors || {A~^—AJ^^Trh)f\\ ■ 

We now apply an exponentially convergent sine quadrature (see, for example, [36]) 
to approximate the integral in ([3]). Since the argument below does not require the 
operator to be discrete, the analysis includes and focuses on the case of 

A-^. 

The change of variable and sine quadrature approximations for A are thus 
A-l^ = r e(i-«//(e/// +A)-i dy 

^ J —oo 
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and 

l=-N 

Here, for any positive integer N, yt := and k := l/y/N. 

To estimate the quadrature error, we start by defining D ^/2 := {z G C, pm(z) | < 
7r/2} and denote by -D,n -/2 its closure. Here uj = arctan(? 7 ) and r] is the index of the 
sesquilinear form H(.,.), see ([9]). For any functions u,v £ and z G C with 

—not in the spectrum of A, we define 

f{z-u, v) := ((e"/ + A)-\, v) . 

Note that fHe(e^) is non negative for z G I?.n -/2 and hence Theorem 12.11 and Re¬ 
mark [2T] imply that f{z,u,v) is well defined for z G D^/ 2 . 

We apply the classical analysis for these types of quadrature approximations 
given in [36] with a particular attention on deriving estimates uniform in u, u G 
L^(n). Using the resolvent estimate fTheorem l2.1|l when iHe(z) > 0 and Remark l2.ll 
when 91e(z) < 0 gives 


{e^I + A)-^\\< 


\e ^|/sin(7r/2 — uj) : for z G Dt^/ 2 , 91e(z) > 0, 

2 _ 

for z G 91e(z) < 0 


Co 


It is a consequence of the above inequality that f{t,u,v) is an analytic function of 
< for < G U^/ 2 - Moreover, for z G Dt^i 2 . 




(42) 


< \\u\\\\v\ 


e-'^^'(^)/sm(7r/2-a;) : for 91e(z) > 0, 

: for 91e(z) < 0. 


Ap(l-/3)«e(2) 


This implies that 

/ CXD 

( 1 /( 2 /- W2;u,u)| -f \ f{y + iTr/2;u,v)\)dy 

-oo 

< (2(^sin(7r/2 - w))"^ -k 4((1 - /3)co)“^) := N{D^/ 2 )\\u\\\\v\\. 
In addition, applying (1^ gives 


j-K/2 

/ \f{t + iy,u,'y')\dy<C, for alH G M. 


Hence, we can apply Theorem 2.20 of [36] to conclude that for fc > 0 

N{D^/2]U,v) ^2 

l=- 


(44) 


nmn,v)dy-k m- 


This leads to the following result for the sine quadrature error. 

Theorem 7.1 (sine Quadrature Error). For N > 0 and k := ^/^/N, let be 

defined by Then for s G [0,1), there exists a constant C{l3) independent of k 
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such that 




(45) 


N{Dt7/2) ■n.2/(4fc) 


2 sinh(7r2/(4fc)) 


ip-/3A 


+ 


+ 




(1 - /3)co 

Similarly, there exists a constant C{(3) independent of k and h such that 


UAf^ - QriAH))nnU^^H^ <Cifd) 


(46) 


N{D^/2) 

2 sinh(7r2/(4fc))' 


= -^V(4fe) 


=,-(l-/3)/fc 


(1 - /3)co 

Proof. Both estimates when s = 0 directly follow from (1441) and the estimates 


“ 1 
oo ^ 

k \f(.^k-,u,v)\<——^eA^-P')/>^\\u\\\\v\\, 

e=N+i ^ k’) 0 

(which are direct consequences of (H^ l. 

For the case s G (0,1), we first consider (|45|) . Using Theorem 14.11 and the 
commutativity of A and (e^/ + A)~^, we have 

^ \\A^/\A-P - Q-/{A))A-^/\ 

Applying (USD shows that the right hand side above equals 

\\A-f>-Q-/{A)\\ 

and so the desired estimate follows again from the s = 0 case. 

For (|46D when s > 0, we apply Lemma [5.21 to see that 

^||(A-^-Q-'^(A,))7r,||. 

The last inequality followed from commutativity and the fact that TTh is a bounded 
operator on iJ®. The result now follows from the s = 0 case. □ 


Remark 7.1 (Critical Case s = 1). The quadrature error estimate still holds when 
s = 1 provided that the discrete and continuous Kato Square Theorems hold. These 
results are given by Remark \f.2\ for the continuous operator A and Theorem 1 6 . 51 for 
the discrete operator Ah. 

Remark 7.2 (Exponential Decay). The error from the three exponentials above 
can essentially he equalized by setting 

(47) Ql^{Ah) := V + 

TT 

i^-M 


TT^/{2k) rv 2pkM Ri (2 - 2P)kN. 


with 
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Thus, given k > 0, we set 


M = 




and N = 


4(l-^)/c2 


1 


1 




and get the estimate 

UA ^- Qk ^ i ^ h )) TTh \\ H>.^Hs < ' 2(1-^)AoJ [sinh(7r2/(4/c)) 

We note that the right hand side above asymptotically behaves like 

1 1 


+e 


l{2k) 


C(/3) 


2/3 2(l-/3)Ao 


,-7rV(2fc) 


0 . 


This, together with the finite element approximation estimates provided in The¬ 
orem 16.21 and Remark 16.21 yields the fully discrete convergence estimate stated 
below. 

Corollary 7.2 (Fully Discrete Convergence Estimate). Suppose that Q and (|8]) 
as well as the elliptic regularity Assumption[J\ hold. Given s S [0,^), set a, := 
+ min(l — 2s, a)) and 7 := max(s -|- a* — /3, s). For d > 7 , then there exists a 
constant C independent of k < 1 and h such that 

Wf G D{A^)nH‘^f 

where Cg.h is given by and is defined by 

In addition, if the continuous and discrete Kato Square Root Theorems hold (see 
Remark \4.2\ and Theorem \6.5\) then the above estimate also holds for s = ^. 

8. Numerical Illustrations for the Convection-Diffusion Problem 


In order to illustrate the performance of the proposed algorithm, we set V = 
Hq{TI) with D = ( 0 , 1)2 and for 6 S R, consider the sesquilinear form 

A{u,v) := / (V u ■ S/ V + b{ux + Uy)v) , for all u, u S V. 

Jn 

This form is regular and the corresponding regularly accretive operator A has do¬ 
main nV. 

In general, it is difficult to analytically compute solutions to u = A~l^ f although 
it is possible in this case. Indeed, we consider the Hermitian form: 


S{u, v) = 


Vu • Vw + ), 


for all u,v gY. 


Fix / G LA{TI) and for > 0, let ic G V solve 

Pl{w, 4>) + A{w, 4>) = (/, 4>), for all 4> gY. 

Putting V = = g-f>(a:-l-y)/ 2 g 0 G V) in the above equation 

and integrating by parts when appropriate, we see that u G V is the solution of 

p{v, 9) + S{v, 9) = 9), for all 9 gY, 

i.e., 

g-&(x+y)/ 2 ^^ -H A)-^f = {(1 + S')-ie-^fo+^)/V- 
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Figure 1. (Left) Decay of e* versus the uniform mesh size ft, in a 
log-log plot. Second order rate of convergence is observed for all 
values of /3. The number of quadrature points is taken large enough 
not to interfere with the spatial discretization error. (Right) Ex¬ 
ponential decay of as a function of the square root of number 
of points M + N + 1 in a semi-log plot (see Remark 17.21 for the 
definition of ft). The spatial discretization is fixed and consists in 
10 uniform refinements of fl. 


Substituting this into ([2]) shows that 

J — ^Hx+v)/2g-!3 ^^-b{x+y)/2 

Since, for example, sin(7ra;) sin(27r2/) is an eigenvector of S with eigenvalue -|- 

672 , 

u := A~^f = -|- 6^/2)“^ sin(7ra;) sin(27ry) 


when 

(48) / = sin(7ra;) sin(27ry). 

The space discretization consists of continuous piecewise bi-linear finite element 
subordinate to successive quadrilateral refinements of D. Figure [T] provides the 
behavior of the errors := \\{A~^ — (,Ah)'Kh)f\\ when the advection coefficient 

is given by 6 = 1 and / is given by (I48p . For a fixed number of quadrature points 

^ h? while for a fixed spatial resolution is exponential decaying. This is in 
agreement with the estimate provided in Corollarv l7.2l 

We now set b = 10, / = 1 and study the boundary layer inherent to convection- 
diffusion problems. For this, we consider 8 successive quadrilateral refinements of 
Q, for the space discretization. In particular, the corresponding mesh size ft := 2“® 
is fine enough for the Galerkin representation not to require any stabilization. The 
value of the approximations A'^^1 for j3 = 0.1,0.3,0.5, 0.7, 0.9 over the segment 
joining the points (0, 0) and (1,1) are plotted in Figure |S] together with the graphs 
of 5 j 4))^1 for P = 0.1 and P = 0.9. The results indicate that the width of the 
boundary layer for convection-diffusion problems remains proportional to the ratio 
diffusion / convection and is therefore independent of p. However, its intensity 
decreases with increasing /?. 
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Figure 2. Approximations of A~^l on a subdivision of the unit 
square using 4^ quadrilaterals and 401 quadrature points. The 
width of the boundary layer for convection-diffusion problems ap¬ 
pears independent of /3 while its intensity decreases with increasing 
p. (Left) Plots over the segment joining the (0,0) and (1,1) for 
B = 0.1, 0.3, 0.5, 0.7, 0.9. (Right) Approximations scaled by a fac¬ 
tor 5 for/5 = 0.9 and/3 = 0.1. 
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